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We consider the evolution of scalar perturbations in a class of non-singular bouncing universes 
obtained with higher-order corrections to the low-energy bosonic string action. We show that pre- 
vious studies have relied on a singular evolution equation for the perturbations. From a simple 
criterium we show that scalar perturbations cannot be described at all times by an homogeneous 
second-order perturbation equation in pre-big bang type universes if we are to regularise the back- 
ground evolution with higher-order curvature and string coupling corrections, and we propose a 
new system of first-order coupled differential equations. Given a bouncing cosmological background 
with inflation driven by the kinetic energy of the dilaton field, we obtain numerically the final power 
spectra generated from the vacuum quantum fluctuations of the metric and the dilaton field during 
inflation. Our result shows that both Bardeen's potential, $(77, fc), and the curvature perturbation 
in the uniform curvature gauge, TL{r\, k), lead to a blue spectral distribution long after the transition. 
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I. INTRODUCTION 



Although current observations 0-0] of anisotropics in the cosmic microwave background radiation (CMBR) strongly 
support the paradigm of scalar field potential-driven inflation as being the source of primordial density fluctuations, 
leading to an almost flat power spectrum 0, [5} , this cosmological model is not short of conceptual problems [|| . 
I One such loophole is the initial singularity which persists in standard inflation Q, but may eventually be addressed 
1 within a more fundamental theory of quantum gravity. String theory for instance admits a class of cosmological 
' models where a period of super-inflation is driven by the kinetic energy of a scalar field. This "pre-big bang" phase 
with growing coupling and growing curvature is expected to last until the background evolution reaches a regime 
of maximal curvature, the bounce, where pure string effects may eventually turn the evolution to an expanding 
Friedmann-Lemaitre (FL) regime. To represent a viable alternative to standard inflation, these cosmological models 
should also reproduce the high-accuracy measurements of the CMBR. To date, the late time spectral distribution 
of adiabatic density perturbations in pre-big bang type universes remains the subject of an intense debate. Some 
authors [8l-lllj recently argued that a scale invariant spectrum (for the ekpyrotic model jjj, iKjf) or a red spectrum 
(for the pre-big bang scenario 0>E1) might emerge from the bouncing regime, while others [l#IT^ favoured a steep 
JL . blue spectrum with n ~ 3 (or n = 4, respectively), ft is thus of prime interest to determine precisely the spectral 
distribution of cosmological models which transit from a collapsing phase to an expanding FL regime. Since this 
discrepancy in the final spectral distribution results from applying different matching procedures to tree-level singular 
1* ■ models, a regularised background evolution could in principle provide the framework to discriminate between these 
. £h ! opposite claims. 

One possibility of smoothing out a curvature singularity of the background evolution is to enhance the low-energy 
effective action of string theory with higher-order curvature corrections (e.g., see ^3-|13|)- The general feature that 
emerges from this scheme is that higher-order corrections generally saturate the growth of spacetime curvature, while 
quantum loop corrections are required to violate the null energy condition and trigger a graceful exit to the post- 
bounce decelerating phase. The inclusion of higher-order corrections in the perturbation equations eventually modify 
the evolution of adiabatic density perturbations 24]. Lately, these perturbation equations have been used |22| to 
study the ekpyrotic model. In this work, the authors have followed the perturbations through a regularised transition 
and have obtained a final n ~ 3 spectrum. Here we shall argue that their method relied on a singular perturbation 
equation. We will then propose a new homogeneous system of coupled linear differential equations involving Bardeen's 
potential <£> and curvature perturbation in the comoving gauge, 1Z. This system remains well-defined at all times. 
Numerical simulations then provide us with post-bounce spectral indices which we find are identical, n$ = n-n ~ 3 
(or n$ = n-n = 4 for the pre-big bang scenario). Thus our result appears in contradiction with the spectral index 



* E-mail: cyril.cartier@physics.unige.ch 



advocated by Durrer and Vernizzi y| and favours the steep blue spectrum put forward in [23j and others. 

The paper is organised as follows. In Sec. [H] we introduce a general action including possible curvature corrections 
up to fourth order in derivatives and we derive the corresponding field equations. Section lllll is devoted to studying 
the classical evolution of adiabatic density perturbations. We first comment on previous results and demonstrate 
that, in our particular class of cosmologies, the perturbation equations cannot be reduced to a set of one decoupled 
homogeneous second order differential equation for each variable. Given a background evolution, we then confront our 
predictions for the late time spectral distributions by numerically integrating a new set of perturbations equations. 
Finally, in Sec. IIVI we recall our main results and discuss possible extensions. For the sake of clarity, we leave the 
details of the calculation to the Appendix. 



II. BACKGROUND EVOLUTION 



As a starting point, we consider the following D-dimcnsional effective action 
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with I a spacetime-length parameter and the higher-order correction terms are given by 

= [ Cl R 2 GB + czGr^v + c^H^^ + c 4 (<P0 ;M ) 2 ] , 
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where R 2 GB = R^\ p ~ 4i?^„ + R 2 is the Gauss-Bonnet combination. Here ^((f)), u>(4>) an d £(</>) are algebraic 
functions of a dimensionless scalar field 0, and with the potential V(<j>) we leave open the possibility of scalar field 
self- interactions. Through the lagrangian density jCS c ', we allow for the inclusion of terms with higher numbers of 
derivatives such as contracted quadratic products of the curvature tensor and define its associated energy-momentum 
tensor by \\f^g T^^Sg^ = 5(\/~g C^°'). Extremising the general effective action Eq. Q yields the covariant 
Euler-Lagrange equations for the cosmological model, 
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where T^ ^ is the energy-momentum tensor of the tree-level terms. The variation of the lagrangian density Eq. 
with respect to the scalar field <f> yields the a' correction A^ which satisfies T^^.^ — Ay<f) ;u according to Bianchi's 
identity. Explicitly, they are given by pH 
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= Cl ^R 2 GB + c 2 G^ [Zrf.^ ~ 2i,^, v - 2£<f>, MV ] 
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Eqs. Jjyl JH| form a complete set of covariant generalised Einstein equations for a background evolution including 
higher-order corrections in the metric and the dilaton field, hence extending their domain of validity to highly-curved 
regimes. For string solutions of physical interest, we wish to restrict to four physical spacetime dimensions, which we 
assume to be described by a FL metric with infinitesimal line element ds 2 — a 2 (77) {—drj 2 + dx 2 ), where r\ denotes the 
conformal time and ' = d/ drj. Inserting this line element in the above covariant equations yields a closed system of 
dynamical equations for the background evolution (cf. Appendix IA 2(1 . Since the quantity encompasses a quadratic 
expression which vanishes in four spacetime dimensions on account of the algebraic identities satisfied by the Riemann 
tensor, we shall neglect it hereafter. 

As an application of these equations, we may consider in particular the pre-big bang scenario of string cosmology 
[T2I Il3l] which requires T — = —uj((f>) and V{4>) — when considered at tree-level in the string frame. In string 
theory, the precise form of the higher-order corrections can be fixed by requiring that our effective action Eq. (JTJ 
reproduces the string theory S'-matrix elements. At the next-to-leading order in the a 1 expansion, this only constrains 
the coefficient of R 2 LV \ p with the result that the pre-factor for the Gauss-Bonnet term has to be ci = 1. But the 
lagrangian can still be shifted by field re-definitions which preserve the on-shell amplitudes, leaving the three remaining 
coefficients satisfying 

£M=e-*, d = l, c 3 = -i[c 2 + 2( Cl + C4 )] . (9) 

The parameter A allows us to move between different string theories : A = 1/4, 1/8 for the bosonic and heterotic 
string respectively, whereas for type II superstrings A = and corrections start at higher order [24I l2fij . Here we shall 
use A = 1/4 to agree with previous studies j2fj- The natural setting C2 = c 3 = leads to the well-known form which 
has given rise to most of the studies on corrections to the low-energy picture. In 0,H(|, the authors demonstrated 
that this set of minimal tree-level a' corrections regularises the singular behaviour of the low-energy pre-big bang 
scenario. Emerging from the asymptotic past vacuum along the low-energy exact pre-big bang solution, the (+) 
branch 0' = (3 + \/3) H with increasing coupling and curvature, the a' corrections drive the evolution to a fixed point 
of bounded curvature with a linearly growing dilaton in cosmic time (the • in the left panel of Fig. ^) . This suggests 
that quantum loop corrections — known to allow a violation of the null energy condition, p + p — would permit 
the crossing of the Einstein bounce 7i e =Tl — (p'/2 |2j| (the subscript -e- denoting a quantity evaluated in the Einstein 
frame) and a graceful exit to a FL decelerated expansion in the post-big bang era, represented by the (— ) branch 
with cf>' = (3 — -\/3) Tt. Little is known about the exact form of quantum loop corrections and we take the freedom 
to parameterise them with an expansion in the string coupling g 2 — of the form £(4>) = + £4 + £,bb^ with 
constant free parameters £a and £b ■ An example of successful exit is shown in Fig. ^ I n general, the combination 
of tree-level a' and quantum loop corrections does not lead to a constant value for the dilaton in a finite amount of 
time. But this can be fixed by introducing a suitable potential or invoking particle production [LsLITi| . 

Given that there exists a class of non-singular cosmologies based on these higher-order corrections, it is then natural 
to investigate the effect of the correction terms on the evolution of primordial scalar fluctuations. 

III. LINEAR PERTURBATIONS 

A. The curvature perturbation 1Z 

To study perturbations, we include small deviations of the background metric, — g^ + Sg^, and for the scalar 
field, cf>. We know from [21j that an homogeneous second-order differential equation for the curvature perturbation 
on uniform energy density hypersurfaces, 1Z, can be derived from the general action Eq. (JJJ including higher-order 
corrections Eq. J5J). To linear order, the decomposition in Fourier modes implies that each comoving wavenumber k 
evolves independently from other comoving modes according to 
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where we used the eigenstates of the Laplace-Beltrami operator, A1Z = —k 2 lZ. Here, we introduced two functions 
of time, z = a^/Q and s. For a minimally coupled massless field cf> and neglecting the contribution of the curvature 




FIG. 1: We illustrate here two non-singular cosmological evolutions relying on the inclusion of higher-order curvature corrections. 

The figure on the left illustrates the flow diagram (<j>,H), where the (cosmic) time coordinate is just a parameter along the 

flow lines. The x-axis corresponds to <f> — (j> — 3H and the y-axis to H = a/a, the Hubble expansion rate in cosmic time in 
the string frame. The (+)/(—) branches correspond to the accelerated/decelerated phases of the pre-big bang scenario without 
potential. The line Eb represents a null Hubble parameter in the Einstein frame, H e = 0. A typical non-singular evolution thus 
emerges from a weak-coupling and low-curvature state, located at the origin of the flow diagram. Then follows a dilaton-driven 
phase with increasing coupling and curvature along the (+) branch and runs into a perturbative singularity, unless higher-order 

corrections force the background to undergo a branch change, cj> — + cj> ^ 0, hence regulating the background evolution. 
Crossing over the Einstein bounce ensures that the Hubble parameter in the Einstein frame is henceforth positive. The figure 
on the right shows the corresponding evolution of H and cf> as a function of the number of e- folds, TV = ln(a). Both regular 
evolutions have been obtained for ci = — f , C2 = C3 = 0, C4 = 1, with a loop control function £(0) = + £a + Cs e< ^- For 
the fast bounce, we set £a = 1 and £b = —2.5 • 10~ 3 . The initial conditions for the simulations have been set with respect to 
the lowest-order analytical solutions at 77 = — 10 6 . A long high curvature regime is easily obtained by trapping the background 
evolution in the a'-asymptotic fixed point (denoted by a •), featuring a constant Hubble rate and linearly growing dilaton 
in cosmic time. This can be achieved by reducing the constant parameters (,a and £s, or by shifting the initial value of the 
dilaton, <j> — > </> — <f>a with, e.g., cf>o = 25. Eventually the string coupling g s = e^^ 2 grows strong enough to allow for a violation 
of the null energy condition p + p ^ and enables a successful exit to the decelerated post- big bang regime. 



corrections, these parameters reduce to Q — > (cf)' /H) 2 and s — > +1. However, in the more general situation where 
the coupling remains unspecified and we allow for higher-order terms of the form of Eq. (J2J), we have (see [2l| for an 
explicit derivation) 
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To be concise, the Qi's and E^s are complicated functions of the background time-dependent quantities, the Qi's 
representing the a' modifications of the tree-level equation. Their explicit form is given in the Appendix. Introducing 
the canonical variable v = zlZ [2^] , the linearised wave equation (|10|l reduces to 
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Through the external potential U (77) , the "pump" field z is responsible for the parametric amplification |29l |3Cj of 
the metric fluctuations. In the case where the background evolution undergoes a period of pole-like inflation (e.g., 
derived from induced grav ity |3l| , scalar-tensor gravity [32] | , the (modified) pre-big-bang scenario 0, IT2I IT3I ] or the 



ekpyrotic model 0,^3H3), the pump field reduces to a power-law, 
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hence U{rj) = q(q — 1)|t/| 2 . Then 



Eq. Ijl3|l corresponds for s = 1 to a Bessel equation and its general solution can be expressed as a superposition of 



Hankel functions of the first and second kind which, in terms of our original variable, reads |2 



fc= : ^£*(*)*tf ) (*M), ^^\l-2q\ . (14) 
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This solution satisfies the usual unitary condition between the coefficients \d,2\ 2 — \di\ 2 — 1, while choosing a pure 
positive frequency state in the asymptotic flat space time for rj — > — oo requires d,2 — 1 and d\ = 0. 

In the large scale limit k/\H\ — > 0, Eq. Ijl4(l reduces to TZ ~ -<4(fc) + ( A;) 1 77 1 1 2<? , where the coefficients A and 
J5 are fixed by the exact solution 114(1. Defining the initial power spectrum by = |a;| 2 A: 3 oc A:™ 1-1 , the spectral 
index carried by each mode before the bounce is thus ua = 3 + 2q and ub = 5 — 2q. For the pre-big bang scenario 
with vanishing potential, we have q = 1/2 and ua — =4. However, adiabatic density fluctuations can also be 
described by the gauge-invariant Bardeen potential in the longitudinal gauge, $> e = & — ^Sip x , which evolves according 
to + 6H$' e + fc 2 $ e = at tree-level in the Einstein frame. In the large scale limit, $ e ~ C(k)\r]\^ < - 1+2qS> + D(k) 
and consistent normalisation yields rip = 1 — 2q = and no = 3 + 2q = 4 for the pre-big bang scenario. Comparing 
the amplitudes of the modes C and D, we clearly have $ e ,C 3> ®e,D right before the bounce, |ry| — > (see [34| for 
the details). Naively, we therefore expect that the dominant spectral index after the bounce remains ric- However, 
since in the radiation dominated era <I> e and TZ differ only by a fc-independent constant, general relativity requires 
them to have the same spectral distribution. Hence the claims that the spectral distribution of TZ is unaltered during 
a bounce and that the dominant pre-bounce spectral index yields the late time post-bounce spectral distribution for 
$ e are contradictory. 

One possibility to discriminate among these incompatible results is to consider a regular background evolution. As 
discussed in plj , the corrections (0 genuinely alter the time-dependence of the variables z = a\[Q and s during the 
high-curvature regime; the source Q of the effective potential for the perturbation may no longer have a monotonic 
growth and may even decrease at the onset or during a high-curvature regime. Apparently, there is also no restriction 
on the sign of the frequency shift occurring when the curvature corrections dominate the background dynamics: s 
may become negative or infinite depending on its particular form. Previous applications of Eqs. (|l(J[) - i|12|) to pre-big 
bang-like scenarios pll have however missed out the singular behaviour of Q during a graceful exit. Indeed the 
components of the Q-denominator E\ — a 2 (T 1 + 2TH. + aQ g ) have to cancel out each other at some particular time, 
since the asymptotic analytical branches of the pre-big bang scenario (T = e - * in the string frame) yields opposite 
signs for Ex : 

-00 cj)' = (3 + V3)H -> E x = -(V3+ l)a 2 He^ < , (15) 

r) s < rj < +00 : tf>' = (3 - VS)H -> E x = (V3 - \)a 2 Ue^ > , (16) 

where rj s represents the time at which higher-order corrections dominate the background dynamics. Furthermore, this 
is the case for any FL spatially-flat bouncing universe which is described by the action Eq. (JIJ, as E\ oc TL — 4>' /2 = 7i e 
has to change sign during the transition between a collapsing phase and an expanding phase. Therefore the pump 
field z and the ratio z 1 / z are singular, which invalidates the use of Eq. (|10|l in the context of bouncing universes of 
the type considered here. 

We now proceed to investigate yet another dynamical equation for the perturbation TZ, for we are interested in 
obtaining the spectrum long after the transition from a regular equation. 

B. Coupled systems (<&, 5<j> x ) and ($,7£) 

The approach we consider henceforth relies on the gauge- invariant Bardeen potentials, $ and ^E", and the pertur- 
bation in the scalar field, 5<p x . Here we shall present our main results, and leave the definition of these variables and 
an explicit derivation of the terms involved to Sec. IA 31 At tree-level, it is well known j^al that these variables satisfy 
homogeneous second-order differential equations at all times. When higher-order corrections are included, however, 
the decoupled wave equation of each of these variables is plagued by a similar singular behaviour during the transition 
regime as the one of TZ = — $ + ^-S(f) x . To see this, we recall that the covariant equations ©-01 yield five redundant 
evolution equations in terms of the gauge-invariant perturbed variables 5" and 5<f> x . They are the four components 
of the Einstein equation, (§), (!•), (■) and the trace-free part ( )tf = (}) — 3 (|) Sj- Those are supplemented by 
the perturbed part of the scalar field equation ((f)). The (°) and (!■) equations are first order in derivatives of the 
perturbed variables, while the (j) and (if) are of second order. Of prime interest is the tracefree part which implies a 
linear background-dependent relation among the perturbation variables 
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where the coefficients c$ and cstj, x are functions of time (see Appendix) , 
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In the limit a' — > 0, we recover 'J + $ = 8<f> x in the string frame, which corresponds to the usual relation ^ e + <f> e = 
when expressed in the Einstein frame. The relation among these variables is modified in a non-trivial way when the 
background dynamics is dominated by the curvature corrections, but the alteration does not exhibit any dependence 
with respect to the comoving wavenumber k. For a regularised background evolution with Eq, ^ at all times (this is 
the case for a wide range of models we have tested), Eq. I|17[) enables us to replace the variable \t = *($, 5<p x ) in the 
perturbation equations. We can then manipulate the remaining equations to extract a set of decoupled dynamical 
equations for the perturbations $ and 5<f> x . However, these equations are plagued by a singular behaviour during the 
transition regime of the kind we previously discussed for 1Z. To demonstrate that we cannot obtain regular decoupled 
second-order wave equations for each of these variables, we focus on the (g) and (°) equations, which completely 
specify the dynamical evolution of the perturbations. These first-order equations can be written as AX' + BX = 
where X = (<£>, S(f> x ) T , and A(rf) and B(r), k) are 2x2 matrices whose coefficients are complicated combinations of the 
background quantities. Provided that the matrix A is invertible, i.e., det(„4) = — 2f2 ^ at all times, we can extract 
an homogeneous system of linear differential equations X' = CX with time-dependent coefficients C(rj, k) = —A~ X B : 
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This homogeneous linear system has a set of n = 2 linearly independent solutions of the form 



X(r ) ,k)=Teiq ? [ / C(f),k)dfj \X( m ,k) 
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where T is a time-ordering operator. If dj (77) are continuous functions and initial conditions are prescribed, the 
system has a unique solution given by Eq. (|20[1 . Although every differential equation of order n can be rewritten as 
a first-order system of n equations, the system viewpoint is more general as we shall recall now. In principle, the 
system Eq. <(T§|) can be reduced to a second-order differential equation by elimination of one of the variable, e.g., 5(j) x . 
Indeed it is straightforward to obtain the homogeneous differential equation 
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But if Cn or C12 are not differentiable, or if C12 = at some particular time 77,, the reduction to Eq. 121(1 will not be 
possible (leaving aside the trivial case for which C12 = at all time) . 

For all background evolutions of the type discussed in Sec.|n] w e have tested, we find that E 2 ^ and ^ during 
the whole evolution, while C12 and C21 change sign. Hence homogeneous second-order equations valid at all times for 
either $ or 8<f) x cannot be derived in our particular class of cosmologies, and we shall use Eq. (|19|l to determine the 
spectral distribution long after the high-curvature regime. Explicitly, the decomposition Cij{rj, k) = Cy (ry) + k Cijirj) 
yields 
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where the coefficients c$(rf) and cs^(rf) are given in Eq. 1|18|) . i s given by Eq. (|llfl and we define 

Sfa) = (^4^10 + 2^5^x1) /(2fi) , 
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Hence, Eqs. I|19|l and 122|l provide us with the necessary tools to determine the spectral distribution of $ and 5(j) x 
long after the high-curvature regime. If numerical simulations were close to infinite precision, we could then deduce the 
evolution of 1Z = — $ + ^p-S4>x an< i ^s late time power spectrum. This in turn would enable us to discriminate among 
the incompatible claims resulting from different matching procedures used in the tree-level analysis. Unfortunately, 
lack of numerical precision inevitably spoils the spectral distribution of 1Z. In a string cosmology background, we recall 
that metric and field perturbations are normalised with respect to the maximal amplified frequency k max = Max(7i), 
for which the perturbation variables have similar amplitudes, \8§\ ~ \^84> x \ ~ Therefore, as long as $ and 8<j> x 

carry red spectra, we will have |<5$| w \&8<t> x I ~ (k m ax/k) 2 \S-iz\ since we expect TZ to have a steep blue k 3 spectrum. But 



for sufficiently small wavenumbers, k <C k max , the spectrum of 1Z obtained from 7Z = — $ + ^7<5</>x will be dominated 
by the imprecision in numerically evolving $ and 6cj> x ; the deduced spectrum of TZ would thus be red, i.e., a measure of 
the inaccuracy of the numerical simulations. To avoid such contamination and erroneous conclusion, we may instead 
consider the variable Y = ($,7?.), whose dynamics also satisfies an homogeneous system of linear differential equations 
of first order, Y' = T>(rj, k)Y : 
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where we have decomposed again the matrix elements according to T>ij (r),k) = T>ij {rj) + k 2 T>ij (77) . Eq. (|25(l is the main 
result of this analysis : It determines the evolution of the curvature perturbation TZ on super-Hubble scales. Indeed, 
the background equations, including those higher-order corrections needed to regularise the background evolution, 
imply at all times the exact cancellation of some coefficients, such that 2?2i = = £>22- Hence TZ remains nearly 
constant on super-Hubble scales, its evolution entering only as a (k/Ti.) 2 correction, while Bardeen's potential may 
evolve drastically on super-Hubble scales. This clearly suggests that, if the perturbation variables TZ and $ are to yield 
the same spectral index long after the transition, the pre-bounce growing mode of the Bardeen potential has to be 
fully converted into its decaying mode during the high-curvature regime. This is confirmed by noting that the leading 
source term of $' on super-Hubble scales, i.e., the coefficient T>u goes from positive to negative during the high- 
curvature regime. The growth of the Bardeen potential is thus turned into a rapid decay during the high-curvature 
regime. 

The complexity of the matrix elements of Eq. I|25|l does not enable us to solve this system of differential equations 
exactly, and we have to resort to numerical integration of the perturbation equations to confirm our prediction about 
the post-bounce spectral distributions. Initial conditions in the asymptotic past can be set according to 
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where we may use the exact tree-level pre-big bang solution Ijl4(l for TZ and its first time derivative. We recall that these 
relations cannot be used at all times since E\ is forced to change sign during the course of the background evolution, 
and we have also shown that there is no second-order homogeneous equation for 1Z valid at all time. However these 
relations are fully adequate to initialise the perturbations $ and 5cj> x at early times rj <C — r/ s < 0. 

Figurc0ilhistrates the results of numerical integration based on the system Eq. (|25H for the pre-big bang scenario of 
string cosmology. We observe that the pre-bounce dominant mode of the Bardeen potential carrying the red spectrum 
is fully converted into the decaying post-bounce mode. Indeed, as shown in [34|. a kink in the spectral distribution 
arises only in the situation where the growing mode of the pre-bounce phase is fully converted into the decaying mode 
after the transition, and one has to wait a sufficiently long time for the decaying mode to decay and the final growing 
mode to dominate. As a result, the pre-bounce subdominant mode yields, long after the bounce, the spectral index 
relevant for the observed anisotropies in the cosmic microwave background; it corresponds to a steep blue spectrum 
with 71$ = 3 + 2q = 4. Hence, although the perturbation variables $ and 1Z carry different spectral distributions 
before the bounce, our numerical simulations confirm that they have the same spectral index sufficiently long after 
the transition when the decaying mode of <3? has died away. This favours the argument put forward in [14j-|l3 1 wne re 
the authors have shown that this is exactly what happens if a matching between the pre- and post-bounce epochs is 
defined by a vanishing jump in the metric and the second fundamental form on the constant energy hypersurface. 



IV. CONCLUSION 



Even though we have not been able to demonstrate from first principles that the pre-bounce growing mode of the 
Bardeen potential $ (which carries a red spectrum) has to be fully converted into the decaying mode during the 
high-curvature regime, we have derived a new master equation which enables us to describe the evolution of adiabatic 
density perturbations on super-Hubble scales. This master equation takes the form of an homogeneous system of 
linear differential equations of first order, explicitly coupling the Bardeen potential and the curvature perturbation 
1Z. On the one hand, we have shown that such system cannot be reduced to second order homogeneous differential 
equations for the class of cosmologies considered here, hence correcting erroneous claims in the literature. On the 
other hand, this master equation clearly indicates that the curvature perturbation is nearly constant on super-Hubble 
scales, its time-dependence entering as a (k/H) 2 correction. This implies that the pre-bounce growing mode of $ 
is fully converted into the decaying mode during the transition. This is confirmed by numerical simulations which 
indicate that the dominant mode of the Bardeen potential carries a blue spectrum long after the transition. As we do 




FIG. 2: We illustrate here the post-bounce spectral distributions of adiabatic density perturbations in the non-singular cos- 
mological background of Fig. We considered the high-curvature transition of short duration in order to reduce numerical 
uncertainty. The picture on the left shows Vx = \x\ 2 k 3 for x £ {TZ, 5<j> x , 5>, $ e } evaluated shortly after the bounce, N = ln(a) ~ 8. 
The wavenumber k is given in units of k max = Max(Ti). The red spectrum attached to the decaying mode of $ is still dom- 
inant on nearly all scales, while TZ carries a blue spectrum. On the panel on the right we compare the spectral distributions 
V x = |a;| 2 fe 3 for x 6 {TZ, $} evaluated at two different times after the bounce : N ~ 8 and N ~ 35. The amplitude of V$ is 
decreasing in time until it reaches its sub-dominant constant mode, which carries a blue spectrum. Hence, both perturbation 
variables $ and TZ lead to the same blue primordial power spectrum long after the transition. Finally, back inside the Hubble 
radius, the spectra start to oscillate and we have $ = ±V37l + est. 



not expect that a smooth transition from the decelerating post-big bang phase to the usual radiation-dominated FL 
epoch (<p = est.) modifies the spectral distribution, we find that the final spectrum of adiabatic density perturbation 
is strongly blue tilted for the class of cosmologies considered here. This result is in complete agreement with j^| 
where the bounce is triggered by a non-local potential instead of higher-order curvature corrections. 
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APPENDIX A 



In this appendix, we provide the necessary tools to reproduce the results discussed in the main text. After intro- 
ducing our notation, we present the different components entering the background equations of motion. Then we 
derive the gauge-invariant expressions which contribute to the linear perturbation equations. Finally, wc introduce 
our matrix notation, and give all the relevant terms in order to derive the systems of differential equations governing 
the evolution of adiabatic density perturbations. 



1. Notation 



To linear order the most general perturbation in a four-dimensional FL spacetime is given by the infinitesimal line 
element 



<4 3) (1 + 2<p) + 2t| 



ds 2 = -a 2 (1 + It) drf - 2a 2 (/% + B t ) dr]dx l + 
Under the gauge transformation 

the induced changes in the perturbed part of the metric read 



; i 2C^|j) -t- 2hij 



dx l dx 3 



(v), 



Ci >—> Ci — 



(?) 



(Al) 



(A2) 



where we use the notation 1 = d/drj and 7i = a' /a. The spatial gauge-invariance of the shear, x = a (P + 7') l— ► X~ a C°i 
enables us to define an infinite number of gauge-invariant quantities. For instance, the Bardeen potentials are built 
exclusively from metric variables 



n 

$ = V — X 
a 



and 



X_ 

a 



Consider a four-scalar (j>(r], x) = 4>{rj) + <5</>(?y, x). Under the transformation laws Eq. (|A1|) . the perturbed part behaves 
as 56 1— > 5c 



b'(P . Hence, one may also construct field-dependent gauge-invariant variables such as 
-7c = <p8<f> = f - —j5(p = --77 5<p v , or 5d x = 5<p-—x = — ~XW ■ 



a a 

as the temporal gauge condition promotes the perturbed variable ip to 



Choosing the uniform-field gauge 54 

become a gauge-invariant quantity, which we denote by tps<t>- Similarly, —^56^ becomes gauge- invariant in the 



uniform-curvature gauge, (p = 0. Up to a sign, it is equal to the curvature perturbation in the comoving gauge 1Z 37]. 
The longitudinal gauge (/? = 7 = 0, hence x — 0) renders the field perturbation 50 gauge-invariant, which we denote 
by S(j) x . We then find 

If we are to consider the longitudinal gauge to fix the temporal degree of freedom, the Bardeen potentials then coincide 
with the perturbations in the metric, $ = ip and \& = r, while <5</> x = <50. 



2. Background evolution 



If we neglect perturbations in the metric and field, the infinitesimal line element reduces to ds 2 — a 2 (j])(—dr] 2 +dx 2 ). 
A solution for the background evolution is then required to satisfy the Einstein equation and the dynamical equation 
for the scalar field, 
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where ^ = d/dcj>, and the different contributions are 
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3. Linear perturbations 



Assuming that there exists a class of background cosmologies which satisfy Eqs. (|A3|) ~ (|A4|) at all times, the equations 
of motion for small perturbations linearised about the background metric are 



6GZ = 

2T 



5T^ u +a'X5T^-2^ l pG^5 ( f> 
Following 35] , the purely geometrical part obeys 
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and so do the terms on the right hand side of Eq. (| A5|) . Using the background solution, the gauge- invariant equations 
of motion for small perturbations linearised about the background metric and field are 
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Here, we defined {T) t f = 7~ l — j^T^S*, i.e., the tracefree part of a tensorial quantity Then, using gauge-invariant 
quantities, the dynamical equation for the scalar field <f> reads 

= 214)5% + 2 (2Huj + wjtf) 5% 
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where we used 



,A-) = (^)) (GI) + (Af))'f. 



The gauge-invariant components entering the above equations are 
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4. Structure of the perturbation equations 

To visualise the structure of these complicated contributions, we write down the perturbation equations in matrix 
form, 



, 



(A6) 



where we use x — (x" , x' , Ax, x) T , x £ {*, *, <50 x }. According to the gauge-invariant contributions of Sec. IA 31 the 
_/Vf O) represent 5x4 matrices, each row of Eq. (| A6|) yielding one of the five component perturbation equations. Here 
we choose that the first three rows correspond to the ($), (!■) and Q) part of the perturbed Einstein equation, while 
the fourth row correspond to the (0) perturbation equations. Finally, the tracefree part of the Einstein equation yields 
the fifth row. Explicitly, after use of the background equations, the matrices entering Eq. I|A6(1 have the following 
structure 
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FIG. 3: Here we illustrate the evolution of the quantities E\, E2 and SI as functions of the number of e-folds, N = ln(a), in the 
context of the high-curvature transition of short duration used to determine the spectral distribution of Fig. [5] As explained 
in the main text, E\ has to change sign in a bouncing universe. The pump field 2 = (</>' E2/ Ef) 1 ^ 2 and z j z are thus singular, 
which invalidates the use of the homogeneous second-order differential equation for the curvature perturbation in the comoving 
gauge, TZ" + 2—72.' + sk 2 1Z — 0. Since E2 and SI do not change sign, the homogeneous systems of linear differential equations 
Eq. 1191 and Eq. 1251 remain well-defined at all times, and can be used to determine the spectral distribution of adiabatic 
density perturbations long after the bounce. 



where we have highlighted the components which are null at all times. It is then obvious that the tracefree part of 
the Einstein equation, i.e., the fifth row of Eq. (|A6|) . implies a linear relation among the perturbation variables 

* = c*$ + C5^5</> x , (A7) 

where the coefficients c$ = —E^E^ 1 and cs(f> = —EgE^ 1 are functions of the background evolution, but do not 
depend explicitly on the wavenumber k. Provided that the background evolution satisfies E 2 7^ at all times, we 
can remove "J from Eq. I|A6() . which yields a system of two second-order differential equations and two constraints 
_/Vf (*)<f> _|_ MS & ^x)Sq> x — 0, the M.^ corresponding to 4 x 4 matrices. It turns out that some of these equations are 
redundant, and one can show that the first two rows, i.e., the (g) and (!•) components of the perturbed Einstein 
equation are sufficient to fully characterise the dynamical evolution of linear perturbations. Explicitly, we get 

f -6E 1 2E 2 c*E 17 \ ( *\ { 2E B 2E 4 -E 10 + c s<Px E 17 
\ E 2 -c*Ex) \ \Ei En-ce^Eh. 

In Fourier space, these first-order equations can then be written as AX' + BX = where X = (<E>, 8<fi x ) T and 

A( ,_f-6E l 2E 5 \ „, ,\_ ( c&E 17 - 2k 2 E 2 c s<t>x E 17 - E w - 2k 2 E A \ 

A[V) -{ E 2 E 4 ) ' - c ^ El E 11 -c 54>x E 1 )■ [AJ) 

Provided that the matrix A is invertible, i.e., det(_4) = —20, ^ at all times, we then have X' — C(r),k)X with 
C = —A^ 1 B. Figure shows explicitly that Q, does not change sign during the evolution for the background evolution 
discussed in the main text. We can then proceed further by replacing the gauge- invariant field perturbation, 8<j) x — 

j£ (<& + TV) , which yields a system of first-order coupled differential equations involving the Bardeen potential $(77, k) 
and the curvature perturbation in the comoving gauge 7?.(?7, k). 

To be complete, we finally present the explicit form of the coefficients of the matrices and M.^ 5 ^ . 

Not all of these quantities are independent and one may use the background equations (including the higher-order 
corrections) for additional simplifications. 
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+c 3 ^ft 2 (2ft' + 3Hft) + 2c^ft 4 ] , 

Q m ee 8a'Xc ia - 2 (,'H , 

Q n ee a'Xa- 4 [-12c 1 ^'H 3 + 9c 2 ^H 2 ft 2 -c 3 ft 3 (2£_'-9£,H)+6c^ft 4 ] , 

Q ee a'Xa' 4 ft 2 [c 2 £ (2H' + H 2 ) + 2c 3 £ {ft' + Hft) + 2c^ft 2 } , 

Q p ee 2a' Xa' 4 [2A Cl £H 3 - 9c 2 £_H 2 ft 2 + c 3 ft 3 - 6^) - 3c 4 £ft 4 ] , 

g 9 ee 2a' Xa' 4 [8ci« (^'W + - ?'W 2 ) - c 2 0' (2^W<^' + 2^W>' + 4^0" - 3£,H 2 ft) 



~c 3 <t>' 2 (ft' + 2£4>" - 2W) + c 4 ^' 4 ] , 



Qr 


= a'Xa 


- 5 [2Ac^'n"H 2 - 3c 2 W0' {i'H<f/ + 2£H</>" + ^H'cj>') 








+c^' 2 + £' {3<j>" - m<f>'} - 6£ {H'tf + 2H(f>"}) - 3c 4 ^' 3 - 




Q s 


= a'Xa' 


6 [24 Cl - n 3 - 9c 2 Z'H 2 ^ 3 + c 3 «/>' 3 (£ V - - 6£'W0') - 




Qt 


= a'Xa 


5 [4 Cl W 2 - - W) - c^' 3 (£' - 3£H) - c 3 /4 (£' - 3£W) 


+ 2c4^' 5 ] , 


Qu 


= a'Xa 


7 [8ciW (e"'W(/)' 2 + ?'</>' {2H'<t>' - H 2 <f>' - 2H<t>"} - £' {H<p'<t>"' + 2H'(f>' 


<j>" -2H4>" 2 -H 2 <j>' >"}) 






-c 2 <p' 3 (2i"U<j } ' + e {2H'4>' + 2H<P" - m 2 /}) - c 3 <j>' 4 (£'>' + W 


- 2W) + c 4 £y 6 ] , 


Qv 


= a'Xa' 


1 [-2A Cl H'H 2 (Z"<t>' - £'<t>") + 3c 2 H<p' 2 (CW + £,"H<f>" + 4£'H'<t>') 





+ C3 0' 2 (-r>' 2 + e'v w - -n + e {xhw + r 2 + w + ew 2 }) 

+3 C4 0' 4 + 3C'0")] , 
Q w = a'Xa' 3 [3c 2 H (£,'H + 2£H') + 6c 4 <f>' + 2£<j>") 

+2c 3 (3£H'<f>' - £<f>" - ?'<t>' + 3W + 3?H<t>')] , 
Q x = a'Xa' 5 [AcxH (2£_"H<j)' - 2£ > 'H<p" + 2?H'<t>' - £,'H 2 cf>') + c 3 <// 3 (-2?<f/ + 3£H(i>' - 2^") 

-c 2 <p' 2 (3?H<t>' + 2£H<t>" + 2£H'<j>' - 3fWV) + 2c 4 £<£' 5 ] • 



